SPECTRAL RECIPROCITY AND MATRIX REPRESENTATIONS OF UNBOUNDED 

OPERATORS 



PALLE E. T. JORGENSEN AND ERIN R J. PEARSE 



Abstract. We study a family of unbounded Hermitian operators in Hilbert space which generalize the 
usual graph-theoretic discrete Laplacian. For an infinite discrete set X, we consider operators acting on 
Hilbert spaces of functions on X, and their representations as infinite matrices; the focus is on e 2 (X), and 
the energy space "Kg. In particular, we prove that these operators are always essentially self-adjoint on 
t 2 {X), but may fail to be essentially self-adjoint on "Hg. In the general case, we examine the von Neumann 
deficiency indices of these operators and explore their relevance in mathematical physics. Finally we study 
the spectra of the fis operators with the use of a new approximation scheme. 
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1. Introduction 

This paper concerns the study of unbounded operators with dense domain in a Hilbert space, and 
their representation in terms of (infinite) matrices. In particular, Theorem 2.6 shows that a "matrix 
Laplacian" is always essentially self-adjoint on ( 2 {X). This class of operators is a generalization of the 
usual discrete Laplacian from graph theory. We also show how the same operator is not essentially self- 
adjoint with respect to the energy space, where the usual t 1 inner product is replaced by a alternative 
inner product given by a natural (quadratic) energy form. We give an axiomatic description of such 
energy spaces H and derive several properties of such spaces from this axiom system. We also prove 
a spectral reciprocity theorem (Theorem 5.16) which establishes an inverse relationship between the 
spectrum of the Laplacian (as an operator on H) and the spectrum of a matrix operator M (as an operator 
on i 2 (X)). The matrix entries of M are defined in terms of a reproducing kernel for H. 

The question of infinite matrix representations of geometric operators arose in a recent project 
[JP09a, JPlOa, JP09d, JP09c, JP09b, JPlOc, JPlOb, JPlOd], where the authors study resistance networks 
and their applications. In these papers, the authors found that that crucial properties of resistance net- 
works may be understood with the use of an associated Laplace operator A, and its various representa- 
tions. The harmonic analysis of resistance forms in the self-adjoint case is worked out in great detail 
in [KigOl, Kig03, Kig09] via an elegant potential-theoretic approach. As noted in [JP09a, Kig03], while 
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tempting to study A as an operator in I 2 , this approach turns out to miss much of the harmonic analysis 
for the given resistance network. The emphasis of the present paper is on the situation where A may 
be only essentially self-adjoint, or possibly even have different self-adjoint extensions. While £ 2 detects 
important spectral data of the Laplacian (and thus also some related combinatorial properties), it sees 
strikingly little of the geometry of the resistance network, in comparison to the spectral theory of A in 
the energy Hilbert space 'Hg, see [JP09d, JP09c, JPlOc]. For a particular problem, the choice of Hilbert 
space {£ 2 or 'Hs or possibly even something else) will play a crucial role in allowing one to extract 
global properties of both the operator and the underlying space (network, graph, or more general set). 
While [JP09a,JP10a, JP09d,JP09c, JP09b, JP10c,JP10b,JP10d] focus on the energy Hilbert space <H 6 , the 
present paper examines the I 2 theory in more depth. Among other things, we examine the deeper reason 
for why A is essentially self-adjoint in t 2 but not in "Kg, and that these two scenarios exhibit drastically 
different boundary conditions in the sense of von Neumann's deficiency indices [vN32, DS88]. 

1.1. Outline. §2 discusses some issues related to the matrix representation of unbounded operators 
with dense domain in a Hilbert space. Special emphasis is placed on a class of operators which we call 
matrix Laplacians, as they generalize the usual discrete Laplace operator on graphs, and can be repre- 
sented in terms of matrix multiplication by an infinite matrix A. Most of §2.2 is devoted to Theorem 2.6 
(and the lemmas required for its proof), in which we prove that a matrix Laplacian A^ acting on l 2 (X) 
is always essentially self-adjoint. 

§3 gives an axiomatic presentation of a class of reproducing kernel Hilbert spaces (RKHS). The 
matrices considered in §2.2 give rise to an RKHS of this type, and the RKHS studied in §4 is a special 
case of this class. 

§4 considers a special case of the reproducing kernel Hilbert spaces of §3 which the authors have 
previously studied in [JP09a, JPlOa, JP09d, JP09c, JP09b], namely, the energy space associated to a 
resistance network. §4.2 describes how A can fail to be essentially self-adjoint (as an operator on "K) 
by explicitly computing an example with deficiency indices (1,1) and giving a formula for the defect 
vector (which is also shown to be bounded and of finite energy). 

§5 returns to the consideration of the matrix M with entries M vy = (v x , v y \ , which first appeared 
as a positive semidefinite function on X x X in §3. The key result in this section is Theorem 5.16, 
which establishes a form of spectral reciprocity between A and D, the diagonalization of M. The exact 
relationship between M and A is made precise in Corollary 5.20; see also Remark 5.25. 

1.2. Basic definitions and facts for unbounded operators on a Hilbert space. In this section, we 
recall some terms and basic results from the theory of unbounded operators on a Hilbert space. This 
material can be found in a standard reference, such as [DS88] or [RS72]. 

Consider an operator T acting on a complex Hilbert space H. We will use T> = dom T to denote the 
domain of the operator T, so T> is always a dense linear subspace of H. 

Definition 1.1. The operator T is Hermitian or {symmetric or formally self-adjoint) iff 

(Tu, v) = (u, Tv), for all u, v e D. 
Definition 1.2. Let T be a densely defined operator in a Hilbert space H. Define 

dom(r) := {v e H I 3C < oo s.t. |<v, Tu)\ < C\\u\\, Vm e dom(r)}. 
In that case, by Riesz's Theorem, there exists a unique w e H such that 

(v, Tu) = (w, u), Vm € dom(T), 
and we set T*v = w. Then T* is the adjoint of T. 

Definition 1.3. If S and T are operators with dense domains dom S c doml c H, then T is an 
extension of S iff T restricted to dom S coincides with S . This is typically denoted S c T, where the 
inclusion refers to the containment of the respective operator graphs. 
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Definition 1.4. An operator T with dense domain D c lis said to be self-adjoint iff T* = T. The 
operator T is said to be essentially self-adjoint iff it has a unique self-adjoint extension. 1 

Definition 1.5. If T is a densely defined operator on H, then T is semibounded iff 



or if the reverse inequality is true. If (1.1) holds, we say that T is a positive semidefinite operator. 

Lemma 1.6. If T is an operator on a Hilbert space, then T is Hermitian iffT c T*. 

Lemma 1.7. Let T be a Hermitian operator on a Hilbert space. Then the essential self-adjointness of 
T is equivalent to 

( i) the closure of T is self-adjoint, or 

(ii) ker(r* ± i) = {0}. 

If T is Hermitian and semibounded, then T is essentially self-adjoint iff 
( Hi) ker(I + 7"*) = {0}, or equivalently, the range ran(I + A) is dense in EL 

Since the Laplace operator T = A discussed below is semibounded, we find it most convenient to use 
criterion (iii). In this case, T is essentially self-adjoint if and only if the eigenvalue problem T*v = —v 
has only the trivial solution v = 0. 

Since the property of semiboundedness is critical in the following, it is shown for the operator A 
acting on the Hilbert space H = ( 2 {X) in Lemma 2.8. The semiboundedness of the operator A acting on 
the reproducing kernel energy Hilbert spaces H = If of §3 is shown in Lemma 3.15. 



We stress the interplay between operators defined on a dense domain in Hilbert space H on the one 
hand and their matrix representation on the other. The questions we address arise only in the case when 
H is infinite dimensional, so we will be considering infinite matrices. Once H is given, we may select an 
orthonormal basis B. Selecting an index set X for B, we note that H is then isometrically isomorphic to 
( 2 {X) = the square summable sequences indexed by X. We will restrict to the case when X is countable, 
i.e., H separable. Our infinite matrices will then have rows and columns indexed by the set X. 

In some of our applications, the set X will be the set of vertices on some weighted graph (G, c) with 
c some (positive and symmetric) function defined on the set of edges in G. In this case, ( 2 {X) will not 
capture the important data for (G, c) and we use a second Hilbert space "Hs defined from an energy form 
for (G, c). In this case, there is a natural Laplace operator A associated with (G, c). It turns out that it 
will have quite different properties depending on whether it is computed in ( 2 (X) or in The matrix 
representations for A will be different for the two Hilbert spaces. Understanding the interrelations of 
these two versions of A in terms of their matrix representations is a main theme of this paper. 

2.1. Matrix representations of operators on ( 2 (X). This paper is primarily concerned with the case 
when X is a countably infinite set, in which case EI = l 2 (X) is separable. Here, ( 2 {X) = ( 2 {X,p) where 
p is counting measure, and we use the usual inner product 



xeX 

and let T : T) — > H be a linear operator on EL For the Hilbert space i 2 (X), we use the orthonormal basis 
(onb) of Dirac masses {S x } XE x given by 



(u, Tu) > 0, for all u e dom(T), 



(1.1) 



2. Unbounded operators on the separable Hilbert space H = t 2 



(X) 





In which case, that unique self-adjoint extension is just the closure of T, in accordance with Lemma 1.7(i). 
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A function uonX will be viewed as a column vector. If A = (a Xd ,) x>ye x is an R- valued function on X x X, 
then Ta(u) = Au is defined by 



Ta : u i-> Au, where (Au)(x) = / j a Xiy u(y) (2.1) 

yex 

(i.e., by matrix multiplication) with the understanding that the summation in the right-hand side of (2.1) 
is absolutely convergent. Henceforth, we describe an object such as A as an infinite matrix with rows 
and columns indexed by X. 

Definition 2.1. The collection of all finitely-supported functions on X is 

cq(X) := {u : X — » C : mL. _ = for some finite subset F c X). (2.2) 

ix\t 

Lemma 2.2. If A — (a X y) x , ye x is an infinite matrix, then matrix multiplication (2.1) defines an operator 



T A : c (X) -» € 2 {X) (2.3) 

if and only if for any fixed y € X, the function x i— > a xy is in t 2 (X). In this case, T A is Hermitian if and 
only ifa Xt y — a y%x for all x,y e X. 

Proof. This is clear because A5 y = a. y is the column in A with index y and (d x , A6 y ) = a xy . The latter 
claim is standard. □ 

Lemma 2.3. Let A — (a XJ ) xyeX be an infinite matrix which defines an operator T A : cq(X) — > ( 2 {X) as 
in Lemma 2.2. Then the following two conditions are equivalent, for two vectors v and w in l 2 (X): 

(i) w(y) — Yjxex a x,yv(x) is absolutely convergent for each y e X, and w e f 2 (X). 

(ii) v e dom(r*) and T* A v = w. 

In particular, tlie action of the operator T* A is given by formula (2.1). 
Proof, (i) => (ii). To show that v e dom T*, note that (Tau, v)p is equal to 



2^ t a u(x)v(x) = 2^jZj a ^y u ^ v ^ = X X a x,y v i x ) u (y) = 2^ w (y) u ( x )> ( 2 - 4 ) 

.vex xeX yeY yeX xeX yeX 

by Fubini-Tonelli. This gives the estimate \(T A u,v)\ < ||w||p||m|^2 by the Cauchy-Schwarz inequality, 
which means v e dom T* . The equality T* A v = w follows from (2.4). 

For the converse, note that w € ( 2 {X) because v e dom(rp. Then the same calculation in reverse 
gives Z x ex u(x)T* A v(x) = Z x ex Zyex a^My)v(x). a 

Corollary 2.4. There exists w € *7Y such that (v,Tau} ( i — {w,u) C 2 holds for all u e dom 7^ if and 
only ifv e domT^ and T* A v — w. If we additionally assume that A is symmetric, the pointwise identity 
(Av)(x) — w(x) holds for all x € X. 

2.2. Matrix Laplace operators on ( 2 (X). In this section, we consider a Laplacian to be the operator 
associated to a matrix satisfying the conditions of Definition 2.5. Our main result in this section is 
Theorem 2.6, which asserts that these three elementary conditions are sufficient to ensure the associated 
operator is essentially self-adjoint, and hence has a well-defined and unique spectral representation. 

Definition 2.5. If X is a countably infinite set, then we say that the infinite matrix A = {a Xty ) xje x defines 
a (matrix) Laplacian iff A satisfies 

(i) 0-x,y = CLyj.i f° r a ^ X 'J 6 

(ii) a Xt y < if x + y; and 

(iii) 2 vex a x,y — 0, for all x e X. 
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In this case, we write for the corresponding Hermitian operator : co(X) — > 1 2 {X) defined by 
matrix multiplication, as in Lemma 2.2. Note that it follows immediately from (ii)-(iii) that a XiX = 
~ ZyexUx} a x,y ^ 0, for each x € X, so the sum in (iii) is automatically absolutely convergent. 

Theorem 2.6 (Essential self-adjointness of matrix Laplacians on l 2 (Xf). If the infinite matrix A = 
(^x,y)x,yex defines a matrix Laplacian on X, then the corresponding Hermitian operator A^ : cq{X) — > 
i 2 (X) is essentially self-adjoint. 

The proof of Theorem 2.6 requires Lemma 2.8, variants of which appear in the literature in different 
contexts, for example, [Kig03, Cor. 6.9] and [FOT94, Thm. 1.3.1]. Theorem 2.6 extends and corrects 
[Jor08, Thm. 3.1] (the result is stated correctly, but there is an error in the proof). 

Remark 2.7. After a first version of this paper was completed, we discovered that Keller and Lenz have 
extended this result to the situation of more general measures in [KL09] and [KL10], as long as the 
measure gives weight oo to infinite paths. (This is true automatically for the counting measure, which 
we use exclusively). Note also that the results of [KL09, KL10] allow for positive potentials (denoted 
therein by c). Consequently, one cannot hope to study the deficiency spaces of A unless one considers (i) 
t 2 spaces with respect to a measure which violates this axiom, or (ii) some other Hilbert space entirely. 
In this paper, we elect to go with the latter option, and hence focus on the energy Hilbert space in §3— §5. 
Related but less general results also appear in [Web09, Woj07]; see also [Woj09]. 

Lemma 2.8 (Semiboundedness of A^ on £ 2 {X)). If the infinite matrix A - (a Xt y) x ,y&c defines a matrix 
Laplacian on a countably infinite set X, then is semibounded and and positive semidefinite with 

(u,Aau)(i = j ^ (— a Xy y)\u(x) — u(y)\ 2 , for all u e cq(X). (2.5) 

x,yeX 

Proof. Note that the right-hand side of (2.5) is a sum of nonnegative terms by Definition 2.5(ii), and 
that it is a finite sum by (2.2). The double summation on the right-hand side of (2.5) is 

2^ a x Ju(x) - u{y)\ 2 = 2j ^ a x Ju(x)\ 2 - 2 ^ a%,y Re(u(x)u(y)) + ^l"^ 2 - 

x,yeX xeX yeX x.yeX yeX xeX 

The last sum on the right side vanishes by Definition 2.5(iii), and similarly the first sum vanishes by 
combining parts (i) and (iii) of the same definition. Thus, the computation continues as 

= -2 2^ a Xi yRe,(u(x)u(y)) = - ^ a x ^u(x)u{y) - ^ a yiX u(x)u(y) = -2(u, A A u)p, 

x,yeX x,yeX x,yeX 

which gives (2.5). In view of assumption (i), we further get that (u, A^w)^ > for all u e co(X). Hence, 
the operator A^ is semibounded and positive semidefinite. □ 

Definition 2.9. An exhaustion of X is a sequence of finite subsets {Ff,}™ =l satisfying F ^ c f t+1 and 

We now return to the proof of Theorem 2.6. 
P roof of Theorem 2.6. Assume that some v e ( 2 {X) satisfies 

2j a x , y v(y) = -v(x). (2.6) 

yeX 

By applying Lemma 2.3 and Lemma 1 .7, we must prove that v = to complete the proof of Theorem 2.6. 
First, observe that (i)-(iii) imply that each of the following functions on X x X is summable, i.e., is in 
e\XxX): 

ax, y \v(x)\ 2 , a x Jv(y)\ 2 , v(x)a Xiy v(y), and a X y\v(x) - v(y)\ 2 . 
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Note that with (i)-(iii), Fubini's theorem applies to the double summations of each of these functions. 
Pick an exhaustion {Fk\T =1 as in Definition 2.9, and then (2.6) gives 



lim 



v + Z a *o> v 00 

yeF t 



= 0, and lim Y a XJ v(y) = -v(x), Vx e X. (2.7) 
The argument in the proof of Lemma 2.8 now yields the following: 



Z X/- fl ^ v W - v( yrf = 2 Z Z ^K^OO - J] Iv(x)| 2 ^ a v - 2 |v(y)| 2 J] a XJ . (2.8) 

.veZjGFt xeXyeF k xeX yeF k yeF k xeX 

Combining (iii) with (2.7) and Fatou's lemma, we can pass to the limit in (2.8). To compute this limit, 
note that for the first term on the right-hand side in (2.8), equation (2.7) gives 

2 2 vOOa^vOO = Z vW Z a *» v W " Z |v(x)|2 = " l|v|1 ^ (2 ' 9) 

.veX yeFt aeX yeFt xeX 

The second term on the right-hand side in (2.8) vanishes because lim^ M 2yeF t fl A,r = 0, by (iii). Con- 
sequently, one obtains the identity 

Z Z ( - fl -«' )|v(x) - v( y^ 2 = - 2||v|1 ?- (2 - 10) 

aeX yeX 

Since the left-hand side in (2.10) is nonnegative (as noted initially) and the right-hand side is nonposi- 
tive, it must be the case that \\v\\^ = 0, whence v = 0. □ 

For future use, we note the following corollary which follows easily from a known characterization 
of positive semidefinite infinite matrices. 

Corollary 2.10. Suppose the infinite matrix A — (a X j) xje x defines a matrix Laplacian on X. If{FkY£ =x 
is an exhaustion ofX as in Definition 2.9, and A(Fk) :— (a x ,y)x,yeF k is the finite submatrix of A corre- 
sponding to Ffo then det A(F > for every k. 

3. Axioms for a reproducing kernel energy space 

In this section, we give some axioms for a certain type of reproducing kernel Hilbert space that distill 
the essential properties of the energy space "Hs discussed in §4. 

3.1. The axioms. Let us fix a set X and suppose that we have a quadratic form Q defined for functions 
u on X with domain dom<3 = \u I Q(u) < oo). Suppose that 'H = dom(2/ker<2 is a Hilbert space with 
respect to the inner product obtained from Q by polarization, that is, under 

(u, v)<h := \ (Q(u + v) - Q(u - v) + iQ(u + iv) - iQ(u - iv)J , (3.1) 
and that H satisfies the following axioms. 
Axiom 1. The constant function \(x) = 1 is an element ofkerQ. 

Axiom 2. For each x e X, the Dirac (point) mass 5 X is contained in dom<2, where 5 X is defined by 

y = X ' (3-2) 
10, y + x. 

Axiom 3. For every pair of points x,y e X, there is a constant C — C xy such that 

\f(x) - f(y)\ < CH/llw, for all f e domQ. (3.3) 
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Remark 3.1. In most cases, it will not be necessary to distinguish between a function u e dom (3 and its 
corresponding equivalence class in H. However, whenever it is useful to make the distinction, we use 
the notation [u]q to indicate the equivalence containing the function u defined on X. 

Remark 3.2. The axiom system above is very similar to the notion of resistance form as developed 
in [KigOl, Kig03] (see also the references therein), although the axioms above evolved independently, 
and from different considerations. Axioms 1-3 allow for slightly more generality than resistance forms. 

Definition 3.3. For a Hilbert space 'H of functions on X, a reproducing kernel is a family {v x } x€ x £ < H 
satisfying 

(v x , u)<h - u(x), for all x e X and for any u e dom<3, (3.4) 

and a relative reproducing kernel is a family {v X y} xje x £ 'H satisfying 

(v Xi y,u)<n = u(x) - u(y), for all x, y e X and for any u e dom Q. (3.5) 

Lemma 3.4. Axiom 3 ensures the existence of a relative reproducing kernel for < H. 

Proof. First, note that (v. c>v , u)^ means {v xy , [u]q)^, as in Remark 3.1. Next, Axiom 3 asserts continuity 
of the linear functional L xiy : ti — > C defined by L x , y u = u(x) - u(y), so Riesz's lemma gives a v x>y e < H 
satisfying (3.5), for each x,y e X. a 

Henceforth, it will be convenient to fix a reference point o e X to act as an origin and consider the 
singly-indexed family {v x } xe x Q ti, where v x = v A . All results will be independent of the choice of o. 

Definition 3.5. Define the possibly unbounded (abstract) Laplace operator with domain 

dom A := span{l, {v x } xeX \u>}} £ < H, (3.6) 

by the pointwise equation 

(Aw)(*):=<*„w>*. (3.7) 

Remark 3.6. In (3.7), the notation (d x , w), H really means <[5 v ]q, w), h , but we can suppress the equivalence 
class notation because any two representatives differ by an element of ker<3. 

Corollary 3.7. dom A is dense in 'H. 

Proof. Suppose that (v x , u). H = for all x e X. Then by (3.5), u must be constant. □ 

Remark 3.8. It is often the case that Aw = 5 X does not have a solution in <H (this is explored in [JP09a]. 
However, Aw = S x - 6„ always has a solution; this follows from Lemma 3.9, just below, and is due in 
some sense to the "balanced" nature of 5 X - 6 y ; see [Soa94, §111.3]. For either Aw = 5 X or Aw = 5 X - 5 , 
the solution w is nonunique precisely when ker A n < H is nontrivial. 

Lemma 3.9. For each x + o, one has Av x (y) — (S x — 6 )(y), for all y e X. 

Proof. From (3.7), we have Av x (y) = ^[<5 v ]<a, v x \ , where [6 y ]Q e "H is the class containing the function 
by defined as in (3.2); see also Remark 3.1 and Remark 3.6. The result now follows via (3.5) by 

{[Syte, v,), w = S y (x) - 5 y (o) = 5 x (y) - 5 (y). □ 

Remark 3.10. From Lemma 3.9, Axiom 2 implies that Au e dom<3 and hence Am represents a unique 
element of Iri. Thus, expressions like (u, Av) w are well-defined, and in particular, so is (u, Av x ), H for any 
x e X, if u g 'H or u e dom<3. 

The following lemma was suggested by (and due to) the referee, for its use in Lemma 3.21. 

Lemma 3.11. Under Axioms 1-3, the set {v x } xe x\{ } is linearly independent. 



8 



P. JORGENSEN AND E. PEARSE 



Proof. Suppose one has a linear combination u = 2£t v .v = where at most finitely many of the 
coefficients % x are nonzero. Then u e domA and Lemma 3.9 gives 

= Au = J] f,Av» = J] 6(5, - 5 C ), 

whence 6 = for all jt e X \ {o}. a 

3.2. Some basic properties of the abstract Laplacian. In this section, we show that the definitions 
given above are sufficient to prove that A is Hermitian and even semibounded. Throughout this section, 
we abuse notation as described in Remark 3.6 and denote both a function and the equivalence class 
containing it by the same symbol. 

Lemma 3.12. If6 xy is the Kronecker delta, then 

(v„ Av v ) h = 6 xy + 1 - 5 xo - 6 yo , Vx, yeX. (3.8) 

Proof. Note that Av y e *H by Remark 3.10, and so 

(v*, Av v ) w = (Av v )(x) - (Av y )(o) = (s x , v y ) H - (s o , v v ) h , 
by (3.5) and (3.7). Again using (3.5), the result follows via 

(d x , v v ) w - (d , v v ) w = (5*00 - 5 x (o)) - (5 (y) - 5 c (o)) = 6 xy + 1 - 6 xo - 6 yo . □ 
Lemma 3.13. The operator A is Hermitian on H. 
Proof. Note that (3.8) is symmetric in x and y, and R-valued. Thus 



□ 



(AV X , V y ) w = (V V , AV A ) H = 6y X + 1 - 5 XO - 6y = 5 X y + 1 - 6 XO - Syo = (V X , AVyj^ . 

Lemma 3.14. The action of A on dom(3 passes to the quotient: [Au]q — A[u]Qfor any u e dom<3. 

Proof. This is equivalent to showing that the kernel of Q is invariant under the action of A. Suppose 
that ifr = YizeF%z v z is an element of ker<3, and that F is finite. Then (ifr,<p)^ = for every <p e fi, so 
with ip = Av x (which is well-defined by Remark 3.10), Lemma 3.13 gives = (t//, Av x ), H = (Aif/,v x ), H , 
for every x e X, so that Aif/ e ker<3 by Corollary 3.7. The conclusion follows. □ 

Lemma 3.15. The operator A given in Definition 3.5 is semibounded as in Definition 1.5. 

Proof. Hue dom A, then u = 2 a£ f % x v x for some finite F QX\\o}by (3.6) and 



< M ,A«).„ = ^y( Vx ,AVy) H = £v6(5,V + 1) = Yj l^l 2 + 

x,y€F x,yeF xeF 

by Lemma 3.12. 



> 0, (3.9) 



Remark 3.16. In fact, one can draw a much stronger conclusion than just semiboundedness from 
Lemma 3. 15: note from the proof that (u, Au), H = implies 2 |£ t | 2 = and thus u — 0. 

Lemma 3.17. Fix y e X and consider if(x) := (v x ,v y j^ as a function of x on X. Let A x denote the 
application of A with respect to the x variable. Then 

A x (v*, v y ) H = (Av x , v y ) M + 6 xo - 1 . (3.10) 
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Proof. Note that tp(x) = v y (x) - v v (o) for each fixed y, so that <p = v y in fi. Then ^ e dom A and 

K (v x , V y ) w = A a .(v v (jc) - v v (o)) = A x v y (x) = 5 y (x) - S y (o) = 6 xy - 8 yo . 

Now (3.10) follows by Lemma 3.12 and Lemma 3.13. □ 

The authors are grateful to the referee for suggesting the above streamlined version of the proof. 

3.3. Foundations of reproducing kernel Hilbert spaces. This subsection aims to give some brief 
historical context for §3 in general, and Lemma 3.21 in particular. 

Definition 3.18. One says M : X xX — > C is a positive semidefinite (psd) function iff 

X TxM(x,y)t y > 0, V{ = tf x } x€X , (3.11) 

xeF 

whenever F c X is, finite. Informally, we describe this condition by saying "M is psd on X". Similarly, 
one says M : X x X — > C is positive definite (pd) iff the inequality in (3.11) is strict for all finitely 
supported nonzero sequences c. 

The theory of positive (semi)definite functions is broad and powerful (see, e.g. [BCR84]) but we 
are interested primarily in two closely related theorems stemming from the work of von Neumann 
and Kolmogorov. The first one (Theorem 3.19) is a generalization and amalgamation of some results 
of [JP09b, §5-6]. The second one (Theorem 3.20) adds the slightly stronger hypothesis of pd (instead 
of psd) and is able to draw a much stronger conclusion: one is able to produce a Gaussian measure on 
the resulting space. The following result is the foundation for the study of reproducing kernel Hilbert 
spaces as developed by Aronszajn [Aro50] and [PS72]. 

Theorem 3.19. Given a psd function M on X, there exists a Hilbert space H~l with an inner product 
(•,•)« and a function v : X — > H such that 

(i) M(x,y) — (v x , v y J for all x,y e H, and 

( ii) cl span v x = < H. 

Moreover, v : X — > 'H is unique up to unitary equivalence when (i) and (ii) are satisfied. In fact, v x is 
defined to be the equivalence class of M(-, x) under a certain quotient map. 

Sketch of proof. The vector space of all finite linear combinations 2 %xM(-, x) can be made into a pre- 
Hilbert space by defining the sesquilinear form 

(y ax M(-,x),J]b y M(-,y)\ := £ cT x M{x,y)b y , 

VreF yeF ' x.yeF 

where F is a finite subset of X containing the support of a and b. One can verify that this satisfies a 
generalized Cauchy-Schwarz inequality, and one can therefore obtain a Hilbert space by modding out 
by the kernel of M and taking the completion. □ 

Theorem 3.20 is an alternative approach to this construction (see [PS72]) which allows one to realize 
the Hilbert space 'H of Theorem 3. 19 as L?{Q., P). This version is more probabilistic in flavour; in fact, 
Kolmogorov's consistency construction is lurking in the background. 

Theorem 3.20. Given a psd function M on XxX, there exists a probability space (Q, P) and a collection 
of random variables [X x ) xe x such that for all x, y e X, 



E(^) = and E(X x ,X y ) = M(x,y). 
Moreover, if M is pd, then P can be taken to be Gaussian. 



(3.12) 
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Lemma 3.21 can be considered as a (somewhat trivial) converse of Aronszajn's theorem, and will be 
useful in §5. 

Lemma 3.21. Given any function v : X — > H mapping X into a Hilbert space, the function defined by 
M(x,y) := (v(x), v(y)) w is pd on X. 

Proof. If £ = {^x\xex is n °t identically 0, then for any finite F c X, 



where w e "H is the function defined by w = Xxef fxV(jf). Note that the final inequality is strict by 



In this section, we introduce the setting of a resistance network (G, c). There are a couple of different 
(but very natural) Hilbert spaces of functions defined on such a domain, both of which are important 
for understanding the underlying network. The study of a network is inextricably linked to the study of 
the associated Laplace operator: note that if A is the adjacency matrix of a network, then as matrices, 
A = cl - A; see (4.6). 

This section aims to compare the ( 2 (G) theory of A (as discussed in §2) with the behavior of A 
on a second Hilbert space of functions naturally associated to the network: the energy space < H&, see 
[JP09a, JPlOa, JP09d, JP09c] and also the references [KigOl, Kig03, LP10, Soa94]. 2 It is defined in 
Lemma 4.4 from an energy form S on functions on (G, c) defined in Definition 4.2. 

The results of §2 imply that the network Laplacian is essentially self-adjoint as an operator on ( 2 {G), 
i.e., on l 2 {G,p) where p is counting measure; see also [KL09, KL10]. However, the action of the 
Laplacian on *Hg is markedly different. In particular, it is not always essentially self-adjoint as an 
operator on < Hs, in sharp contrast to Theorem 2.6. Example 4.8 illustrates this phenomenon with an 
explicitly computed defect eigenvector and (nonzero) deficiency indices. 

It also turns out that there is no natural onb for fisl the natural candidate would be the Dirac masses 
{8 x }xeg, but these are not orthogonal and typically don't even have dense span in 'Kg. Consequently, we 
rely on a reproducing kernel \v x } x£ g\{o), as developed axiomatically in the previous section. In fact, this 
is part of the motivation behind §3. 

Due in part to their close relation with Markov chains, there is a massive literature on resistance net- 
works (not always using this terminology). Many studies use Hilbert space techniques, but almost all of 
these focus on i 2 (G,p); see [Soa94,Chu01] and the references therein; other articles use methods from 
potential theory and discrete harmonic analysis [KigOl, Kig03]. See also [LP10, §9] for an alternative 
view on the energy space, presented in terms of an t 2 space of functions on the edges of G. 

4. 1. Networks and the energy space. 

Definition 4.1. A resistance network is a connected weighted graph (G,c). Here G = (G°,G') is a 
graph with a countable vertex set G°, and at most one edge e e G 1 between any two vertices. From 
this point onward, we write x e G to indicate that x e G°. The adjacency relation on G is determined 
entirely by the conductance function c : G° X G° — > [0, oo), a nonnegative and symmetric real-valued 
function denoted c„. = c(x,y). We say x,y e G are connected by an edge of weight c„ if and only if 
c xy > 0; in this case, we write x ~ y. Vertices may not have finite valency, but they must have finite total 
conductance: 



(3.13) 




Lemma 3.11. 



□ 



4. The Laplacian as an operator on the energy space 




(4.1) 



•Hg is different from the space D discussed in [Soa94] (also called in [Kig03, Prop. 2.19]), but the two are closely 

related; see [JP09a, §4.1], for example. 
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We also assume c xx = for every x e G. 

In Definition 4.1, the term connected means that for all x,y e G, there is a finite sequence {x = 
xq, X\, . . . , x„ = y} c X such that c XiXi _, > for i = \,...,n. There is a bijective correspondence 
between the class of resistance networks and the class of irreducible reversible Markov chains; the 
correspondence is given by considering G° as the state space and defining the transition probability by 
p(x,y) = c xy /c(x), for vertices (states) x and y. 

Definition 4.2. For functions «,vona resistance network, one can define the (sesquilinear) energy form 

8(u, v) = i J] c -n'(«W " «00)(v(*) " v(y)) (4.2) 

x,yeC 

with domain domfi := \u : G — » C I £(w, u) < oo}. One says that S(u) :- S(u, u) is the energy of m. 

It is clear from (4.2) and the connectedness of the network that &(u) — iff u is constant, so kerfi = 
CI. Therefore, we define an equivalence relation by u ~ v iff u(x) - v(x) = k for some fixed k e C. 

Lemma 4.3. Under the above equivalence relation, and with || • He = V£(s 0» tne quotient 

<K e := = {w + CI ! w : G — » C and \\u\\ s < oo) (4.3) 

kero 

is a Hilbert space, and the elements ofH& are functions on G modulo constants. 

Proof. It can be checked directly that the above collection of (equivalence classes of) functions on G is 
complete via an isometric embedding into a larger Hilbert space as in [LP 10, MYY94] or by a standard 
Fatou's lemma argument as in [Soa94]. □ 

Definition 4.4. The energy space is the Hilbert space with inner product (u, v)g := &(u, v). 

Theorem 4.5. The energy space is a special case of the axiomatic presentation in §3. 

Proof. Note that (4.2) gives 



(s x ,s y ) e 



and 



(6 X , S x )s = c(x). 



(4.4) 



In particular, the condition c(x) < oo ensures i5 A e *Hg for every x e G, and so Axiom 2 is satisfied. To 
see that Axiom 3 is satisfied, one can argue as in [JP09a, Lem. 2.4]: since G is connected, choose a path 
{*;}"_„ with xq = y, x„ = x and c XuX ._, > for i - and the Schwarz inequality yields 



\L x , y u\ = \u(x) - u(y)\ 



Z 



(u(Xi) - u{Xi-i)) 



< k 2 £(u), for k ■■ 



1/2 



Consequently Lemma 3.4 applies and we have a relative reproducing kernel {v x } x ec e ^fe> as m Defi- 
nition 3.3, given by v x := v x>0 . Although the elements of fis are equivalence classes, computations can 
be performed using representatives whenever these computations are independent of the choice of rep- 
resentative. Abusing notation, we may take the function u to be the representative ofae 7Yg satisfying 
u(o) = 0. 3 □ 



; After an initial draft of this paper was complete, we discovered that researchers studying metrized graphs use a similar object; 
in [BR07,BF06] this is called the "J-function" and is roughly given by j z (x, y) = v y z (x). The two objects do not precisely coincide 
because for metrized graphs, x,y,z may be points in the interior of a edge, as edges are isometric to intervals in that context. 
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Remark 4.6. Since one may add a constant function without changing the energy, dom<5 = "H ffi C. 
Then, as in [LP10, Ex. 9.6(b)], one has 

= (4-5) 

C 

Upon combining Definition 3.5 with (4.2), one obtains the Laplacian as the (graph) closure of the 
operator defined pointwise on the dense domain dom A = spanjl, {v^} xe y\{ }} by 

(Au)(x) = c xy (u(x) - u(y)). (4.6) 

y~x 

Remark 4.7 (The meaning of Au). Note that Au is a function on G, not an equivalence class of functions 
(the differences in (4.6) specify the value of Au{x) unambiguously). 

It is also clear that A is Hermitian on fig; note that Corollary 3.13 holds in this context. It is also the 
case that A commutes with conjugation, and this ensures that the deficiency indices of A on "He will be 
equal. §4.2 discusses a situation in which A on 'Kg has deficiency indices (1,1). 

Using the standard onb {5 x } xe c for £ 2 (G), and the matrix A with entries a c v = —c xy , formula (4.6) is 
equivalent to matrix multiplication: 



Au = Au, (4.7) 

so that a x y defines a matrix Laplacian on ( 2 {X) in the sense of Definition 2.5. In fact, the only real 
difference between Definition 4.1 and Definition 2.5 is the addition of the connectedness condition, 
which appears in this section to ensure that the kernel of the energy form contains only (globally) 
constant functions. 

4.2. The Laplacian can fail to be essentially self-adjoint on < H £ . 

Example 4.8 (The geometric integers). For b > 1, consider the network (Z + ,b n ) consisting of the 
nonnegative integers with an edge of conductance b k connecting the vertex k — 1 to the vertex k: 



12 3 

See [JP09c, JP09a, JP1 0a, JP09d, JP09b] . 

Proposition 4.9 (Defect on the geometric integers). As an operator on the energy space of the network 
(Z + , b"), the Laplacian is not essentially-self-adjoint. 

Proof. We will explicitly construct a function u which has finite energy and which satisfies Au(n) = 
-u(n) at every vertex n in the network. To this end, recursively define a system of polynomials {ip„} and 
{if/ n ) in the variable r by 



1 



1 

1 +r" 



1 



1 

1 +r 2 



1 

1 +r 



(4.8) 



Putting r = I, the desired function u is defined by u(n) := i//„(l/b). Note that ip„ - ip n -\ + ii n -\ and 
*A« = ifrn-i + i^fn- Hence 

u(n) - u(n - 1) = i/r n (l) - ^„_!(i) = «A„-i(i) + r'V„(i) - 1^-1(5) = ^n(|) 
and therefore, suppressing the evaluation at the fixed value r = l/b,we have 

tp n = b n (u(n) - u(n - 1)) (4.9) 

Consequently, Au(n) = (p„ - tp n+ \ = -if/„ = -u(n) implies that Au = -u. The proof will be complete 
once we show that u e 'Kg, which is carried out in Lemma 4.1 1. □ 
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We will need the following lemma for the proof of Lemma 4.11. 
Lemma 4.10. There is an m such that 

tp„<n m , and tf/„ < (n + 1)'" - n m forallneZ + , (4.10) 
where ip„ and tff n are the polynomials defined in (4.8). 

Proof. We prove both bounds simultaneously by induction, so assume both bounds of (4.10) hold for 
n - 1. The estimate for ip n = ip n -\ + ty n -\ is immediate from the inductive hypotheses. For the tff n 
estimate, choose an integer m so that 



m(m - 1) > maxffV • t > 0} = 
Then n 2 r" < m(m - 1) for all n, so 



2 ^ 2 



elogb 



m(m -1) In - 1 \ In + 1 
2 + r" < 2 + — — r — - < + 1 



by using the binomial theorem to expand ( — I = [1 ± - I . Multiplying by n m gives 

(n m - (n - l) m ) + r"n m < (n + 1)"' - n m , 
which is sufficient because the left side is an upper bound for t// n = i/f„_i + r"yj„. □ 

Lemma 4.11. The defect vector u(ri) :— ij/ n (^) has finite energy and is bounded. 

Proof. Applying Lemma 4.10 to the formula for £ yields 

oo oo oo 

&(u) = J] b"(u(n) - u(n - l)) 2 = J] r >« ^ Z r "" 2m = U - 2m(r) < £X '' 

n= 1 n=l n=l 

since a polylogarithm indexed by a negative integer is continuous on R, except for a single pole at 1 (but 
recall that re (0,1)). 

To see that u is bounded, combine (4.9) and (4. 10) to obtain b" (u(n) - u(n - 1)) < n m , for some fixed 
m, whence the sequence of increments is summable in much the same way: 

oo oo 

lim u(n) - u(0) = Y (u(n) - u(n - 1)) < ^ r n n m < oo. □ 

Lemma 4.1 1 ensures that the defect vector is bounded; in the example in Figure 1, the defect vector 
has a limiting value of * 4.04468281, although the function value does not exceed 4 until x = 10. The 
first few values of the function are 



_ [3 17 173 3237 1 14325 7774837 1032268341 270040381877 140010315667637 I 
~ [2' 8 ' 64 ' 1024' 32768 ' 2097152' 268435456 ' 68719476736 ' 35184372088832 ' ■ ■ • J 

* [1.5, 2.125, 2.7031, 3.1611, 3.4889, 3.7073, 3.8455, 3.9296, 3.9793,4.0080, . . . ] 
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4.044 - 
4.043 - 
4.042 - 
4.041 - 
4.049 - 



FIGURE 1 . A Mathematica plot of the defect vector u of Example 4.8 and Lemma 4.1 1, The left plot 
shows u(x) for x = 0, 1, . . ., 10, and the plot on the right shows data points for u(x), x = 10, 1 1, 12, 



5. Finite approximates 

As mentioned in the previous section, when considering A as an operator on a reproducing kernel 
Hilbert space fi, it is not possible to use the matrix representation of (4.7) because {6 x } x€ x is not an onb 
for <H. Therefore, we change to a different representation of H as laid out in §3. 

In this section, we return to the setting of §3, where X is any (infinite) set, Q is a quadratic form on 
functions on X, and 'H = dom<2/ker<3 is a Hilbert space with (relative) reproducing kernel [v x ] x€ x. For 
studying infinite sets X, it will be helpful to consider a filtration by finite subsets, partially ordered by 
inclusion. With this aim, we pick a finite subset F c X and study the "restriction" of M and functions 
u to this subset. Note that we do not restrict the support of the functions under consideration: we 
restrict the index set of the representing functions {v x } xe x, in the spirit of Karhunen-Loeve; see [Ash65]. 
This is akin to using cutoff functions as Fourier multipliers, and leads to a form of spectral reciprocity 
between the associated Laplace operator, and its "inverse" M, in the sense described in §5.1. The exact 
relationship between M (actually, its diagonalization D) and A is made precise in Corollary 5.20; see 
also Remark 5.25. The application we have in mind is a resistance network as discussed in §4 but all 
results are phrased in the context of the reproducing kernel Hilbert space of §3 so as to keep the scope 
of discussion more general. 

Definition 5.1. Let *V := span{v A } A(E x\{o) and *V(F) := span{v A I x e F}. We also write ( 2 (F) for the 
subspace of functions in ( 2 {X) whose support is contained in F. This may seem trivial when F is finite, 
but the notation helps distinguish between the two different inner products in use. 

Definition 5.2. Define <3> : t 2 {X) -> "H on dom$ = span{<5 A } AE jr by <M A = v x . 

Remark 5.3. The operator <t> is typically not closable. To see this, we show why the adjoint is not 
typically densely defined. First, pick £ e span{(5 A ) and u e <V, and compute <$>*: 

<f , <s>*u) e = U )<h = J] ^< v - m >« = Z t~ x(ll(x) ~ u{o)) - 

xeX xeX 

So for an equivalence class u e ( H, note that <b*u is the representative of u that vanishes at o. For u e "H, 
let us denote by m (0) the representative of u specified by u(o) = 0, so that <E>*m = m (0) . Then 

dom0>* ={ue<H\ « (0) e £ 2 (X)}. 

It is easy to see that this class is not dense in fi; see [JP09a] for examples in the case H = TYg. 

Definition 5.4. For a finite set F c X \ {o}, we have <6£ = Z xeX £(x)v x , for all £ e 'V(F). Define M to 
be the matrix of <f>*<f>, that is, 

M xy = (d x , ®*®6 y ) e = (v A , Vy ) H , Sx,y e X, (5.1) 

and let Mp := M\ Fxf . be the submatrix of M defined by deleting all rows and columns corresponding to 
points x i F, i.e., Mf is an \F\ x \F\ matrix with entries 
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(M F )„ = (v„v,) w , Vx,yeF. (5.2) 
In general, one may have v x e *V(F) with support extending outside of F; examples are given in [JP09b]. 

Note that since £ e ( 2 (F) is finitely supported, 

M&x) = M x ,yt(y) = J] (v„ v,) w f (y) = /v», J] £(y)v y ) = <v„ 0^).„ = 0>£(x) - Q>{(o) (5.3) 

veF veF < yeF ' 

Definition 5.5. Denote the spectrum of M F by A F = {A F .} for some enumeration j — 1,2, ... , |F|. Note 
that A F > by Lemma 3.21 and that M is diagonalizable with eigenfunctions £jj = e £(F ) = ( 2 {F). 
That is, the spectral theorem provides an orthonormal basis (onb) with 

M F % j = X£ j for each ;',F. (5.4) 
For convenience, we often suppress the index and write (5.4) as M F g A = A^ A . 
Definition 5.6. For a finite F c X, and = A£ A as above, define 

«.i := -j= ^ 6W v j- (5-5) 



Lemma 5.7. 77;e operator *¥ F : ^ 2 (F) — > *V(F) defined by *F F (i;j) = is unitary, and consequently 
{u\}xek f is an orthonormal basis in *V(F). 

Proof. For x, y e F, compute 



and since £j is an eigenvector, this continues as 



where 5# is the Kronecker delta, since {£}} is an onb for ( 2 {F). □ 

Definition 5.8. By Lemma 5.7, we may let P F be the projection to span^l^A,,. In Dirac notation, this 
is 

AeA f 

Note that P F is projection to 'V(F). 

Lemma 5.9. With respect to the onb {u A }, one has 

P F v x = ^ A U2 £ A (x)u A , for all x e F. (5.7) 

AeA f 

Proof. Let x e F. Then compute 

,leA f ,leA f /leA F y£F 

by (5.6) followed by (5.5). Continuing, 
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P F v x = V -^-(M F ^ A )(x)u A = ^=&(x)uji, 

AeA F V/t AeAf V/i 

since % A is an eigenvector. Note that A e R + , since M is positive semidefinite by assumption. It remains 
to observe that P F v x = V x for x e F, but this follows from Definition 5.8. □ 

Remark 5.10. In the language of Theorem 3.19, equation (5.7) takes the following form: 

v= Yj V* (5.8) 

where {% A } is an onb for ( 2 (F) and {w^} is an onb for *V{F ). The significance of this symmetric expression 
of v is that it allows us to compute a norm in <H (where the sum would be over x e F ) by instead 
computing an { 2 norm (where the sum is over A € A F ). For an example, see Corollary 5.24. 

In [JP09d], the authors show that for Iri - Iri^ one can construct a Gel'fand triple S B £ "He £ S' B , 
isometrically embed 'Kg e — » L 2 (»Sg,P). Here >S fi is a space of "test functions" which is dense in 7Yg, 
but comes equipped with a strictly finer Frechet topology, and S' s is a space of "distributions" obtained 
by taking the dual with respect to this topology. Elements u e TYg can then be extended to functions 
on S' B via S(£) = (u, £)g for £ e S' E . As P is a probability measure, one can then interpret \v x } xe c as a 
stochastic process, i.e., a system of random variables indexed by the vertices of the underlying graph. 
In this context, (5.8) becomes an instance of the Karhunen-Loeve decomposition (see, e.g. [Ash65]) of 
a stochastic process into its random and deterministic components: 

= Yj ^ ® fi <<^)) " G -^ n (5-9) 

^sA f 

In fact, it turns out that {u A } Ae A F is a system of independent identically distributed Gaussian A^(0, 1) 
random variables, for any finite F c X. See also §5.2 for more relations to Karhunen-Loeve. 

5.1. Spectral reciprocity. In this section, we explore the relationship between M and A. In particular, 
the Spectral Reciprocity Theorem (Theorem 5.16) shows how M and A are (almost) inverse operators, 
and explains why the eigenvalues of M are (almost) the reciprocals of the eigenvalues of A. 

Definition 5.11. Denote the diagonalization of M F by 



(5.10) 



A \F\ 

where P Ul is projection to span{u A } AeAF . Note that D F l is a well-defined operator on i 2 {F) of rank 
\F\ <oo. 

Definition 5.12. Let P^:'H-* < H be the projection of 6 to 'V(F). That is, P° F = \P F 5 ){P F 5 \ in Dirac 
notation. 

Definition 5.13. If {<; A } is the onb of eigenvectors of M F , denote the expectation of % A by 

X £F 

Lemma 5.14. If 5 is a Dirac mass at the origin, the expansion of P F 6 with respect to {u A } is given by 



(5.12) 
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Proof. Using Pp = P* F , (5.5), and the fact that u A e 'V(F), we compute the coefficients: 

(u A ,P F 5„), H = {P F u A ,6 ), H = (u A ,6„), H = V % A (x){v x ,5 ), H = — l — V 

where the last line follows by Lemma 3.9, since x + o. □ 

Definition 5.15. The compression of A to F is the restricted action of the operator A to *V{F ), and it is 
given by PfAPf. 

Theorem 5.16 (Spectral reciprocity). If F c X \ {o] is nonempty and finite, then 

PfAPf = (DD^'O* + P° F . (5.13) 
Proof. For A, k e A F , we have {u A , P F AP F u K ), H = (u A , Au K ), H because u K e 'V(F). Then 

{u A , P F AP F u K ), H = V &W&(y) (v., Av y ) = V WMMtfxy + 1) (5-14) 
by (5.5) and (3.8). The computation of (5.14) continues as 



= ~r= X &(*)&(*) + Z^W-P Z&OO = ~F= Z&W&OO + ^E(£i)^-E(^). 

(5.15) 

Since u A is in domO* automatically for finite F, the right side of (5.13) is 

(h a , (OD^'cD* + = (o'h^D^O'h,)^ + {u A ,P° F u K ), H , 

which matches with the right side of (5.15), by (5.12). This verifies (5.13) on the onb of Lemma 5.7, 
and hence for all of r V(F). □ 

Remark 5.17. We refer to Theorem 5.16 as the Spectral Reciprocity Theorem because it relates the 
eigenvalues of A to the reciprocal eigenvalues of its inverse, on any finite F c X. 

Suppose one writes the matrix for A as in Definition 2.5, so that rows and columns are indexed by 
points of X. Let A be the matrix which results from deleting the row and column corresponding to a 
chosen point o. Corollary 5.20 makes precise the well-known statement that one can invert the Laplacian 
after deleting the row and column corresponding to a point o. 4 In particular, without deleting the row 
and column of o, one is forced to contend with an auxiliary term 1 in (3.8) (which corresponds to the 
projection P" = \o)(o\ to the 1 -dimensional subspace spanned by So). 

Lemma 5.18. For every nested sequence of finite sets {F„} ne m with (J F n = X \ {o}, the limit of P Fn AP Fi 
exists and with dom A as in (3.6), 

A = lim P F AP Fn , (5.16) 
in the strong operator topology, that is, lim,,-^ \\P Fn AP Fn v — Av\\<h for all v e domA. 
Proof. Let / e dom A so that there is some finite set F c X \ {o} for which 



xeF 



Without loss of generality, let {F n }^ =l be an exhaustion of X \ {o} with F c F\. Then P F f - /, and 



^Recall that if M is a Hermitian matrix acting on a finite-dimensional Hilbert space H, then the restriction of M to the 
orthocomplement of the zero eigenspace is invertible. 
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P Fi AP F J = P Fn Af > Af, 

since P Fn increases to the identity operator: 



\\P F AP F J - Afhf = \\(P F „ - DA/|| W > 0. □ 

Definition 5.19. Let P° : = |<5 }(£ | = proj.„[<5 ] be the rank-1 projection on "H defined by {u,P°w). H = 
(u,8 ) n (6 ,w) n . 

Corollary 5.20. The limit A — P° = lim^^, OD^'cD* exists, for any exhaustion {F n } ofX \ {o}. 
Proof. Since arguments exactly analogous to those in Lemma 5.18 give P° p > P°, we have 

lim ftDp 1 ®* = lim P F ,AP Fii + lim Pi. = A - P°, 

by applying Theorem 5.16 and then Lemma 5.18. □ 

5.2. Spectral measures. Recall from Definition 5.13 that E(£y) = Yjxef = £/) > an d tnat from 
Lemma 5.14, the expansion of P F 6 a with respect to {u A } is given by 



V ^ i -n 

Pf5 q = - / j — 7=-**- ( 5 - 17 ) 

AeA F V/l 

Definition 5.21. Since P F AP F = DZ 1 + P F 5„ is the 7 x J matrix T F whose (j, k) th entry is given by 



denote the spectrum of this matrix T F = [t;u\ by S F = spec(7V) = {cr F } J _, . 

Remark 5.22. In Definition 5.21, it is important to note that t and Aj all depend on the choice of 
F. However, for ease of notation we suppress this dependence and also henceforth write crj = a F . 

Recall from Definition 5.5 that A F = spec(Mf) = {Aj}j =1 . 

Corollary 5.23. For any finite subset F c X \ {o} with \F\ = J, one has j Y/j=\ E(£/) 2 = L 

Proof Since X F coincides with the constant vector 1 on F, we use P^u = (^j, uj gj to compute directly 

2 |e(^)| 2 = Z |(*F.6),f = Z l|7 W |2 = M l = \f\ = J. □ 

Corollary 5.24. For any finite subset F QX\{o) with \F\ = J, one has < \\P° F \\ < -jj-j. 

Proof Using Corollary 5.23 and Definition 5. 12, \\P° F \\ = ||P F 5 ||^ = j/ j=x E(£ ; -) 2 = J. See Remark 5. 10. 
Then the double inequality follows by estimating by the largest (but clearly finite) eigenvalue and the 
smallest (but clearly strictly positive) eigenvalue. □ 

Remark 5.25. When A is not essentially self-adjoint, the presence of P° F (as in (5.13), for example) 
makes it impossible to obtain self-adjoint extensions of A via a filtration by finite subsets. This obstacle 
can only be overcome by passing to spectral measures. If dom A is as in (3.6), then the spectral measure 
of some self-adjoint extension of A comes from the weak-* limit of linear combinations of of equally 
weighted Dirac masses: 
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= 7EV '5 IS) 



Here, p F refers to the spectral representation of P F AP F , and we are relying on standard tools from 
the literature. Indeed, approximation of measures with the use of spectral sampling is a versatile and 
powerful tool. For approximation in the weak-* topology on measures (as in the present context), see the 
excellent reference book [Bil99] for details. When applied to spectral measures, these approximations 
were first studied in the book by M. Stone; see [Sto90, Ch. X]. The approach in [Sto90] is especially 
amenable to our present applications: a main theme is the study of unbounded operators in Hilbert 
space, realized concretely as banded infinite matrices. This is illustrated in the following diagram: 

^A^^A (5.19) 



T 



weak-* 



In the limit of (5.18) as F — > X, may jj. F become a smooth measure. The key point is that considering 
the limit of P F AP F as F — > X does not take one far enough. However, consideration of the spectral 
measures p F of P F AP F shows that each weak-* limit ft is the spectral measure of some self-adjoint 
extension A of A, and by general theory, every self-adjoint extension of A arises in this way. 

In the preceding discussion, F — > X refers implicitly to a limit with respect to an exhaustion {F n } n eN< 
where c F^+i and U^Li = X \ {o}. Note that the limit A = lim F ^x Pf^Pf is unique (see 
Lemma 5.18) and hence independent of the choice of exhaustion {F„}„en- However, the nonuniqueness 
of weak-* limits corresponds to the fact that ft - lim F ^ x I->f ma y depend on the choice of exhaustion. 
Different weak-* limits may correspond to different self-adjoint extensions A of A. 

5.3. Spectral reciprocity for balanced functions. Balanced functions are functions which sum to 0. 
In the context of resistance networks (see §4), a balanced function is the divergence of a current flow 
with no transient component; these functions are mentioned briefly in [Soa94, §111.3]. 

Definition 5.26. A function £ : X — > C is balanced iff £ has finite support and Yixex €( x ) - 0- Denote 
the space of such functions by S. For any finite F c X \ {<?}, let < B F denote the collection of functions in 
& whose support is contained in F. 

Recall from Definition 5.1 that *V '■= span{v A } AE x\(o) and 'V(F) := span{v x I x e F], and from 
Definition 5.2 that (D : ( 2 {X) -> <H is given by <5(<5 A ) — v x on domO — spanf^xj^gx- 

Definition 5.27. Denote the subspace of *V with balanced coefficients by 



:= 0(S) = mt) :^S}, (5.20) 

and similarly for 'Vo(F) := <&{B F ) = {<D(£) i £ e S F ). 

The following curious fact can be found in most introductory books on functional analysis. 

Proposition 5.28. Let Abe a topological vector space, and let Aq be a dense linear subspace. Iff is a 
linear functional on Aq, then ker/ is dense in A if and only if f is discontinuous. 

Lemma 5.29. S is dense in { 2 (X) if and only ifX is infinite. 

Proof. Define / : S — > C by /(£) = Yixex&- Note that X is finite if and only if the constant function 
1 is in £ 2 (X), which (by Riesz duality) holds if and only if / is continuous on £ 2 {X). The result now 
follows from Proposition 5.28. □ 

The next lemma indicates how <D "intertwines" the spectral densities of A and M. 
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Lemma 5.30. For all £ e £, one has <<D(£), A<D(£)>« 

«D(0,A<D(0>« 



||0(^)||^, and hence 



mm 



<H 



(5.21) 



Proof. The first identity is immediate for £ e S by (3.9). For the second, note that 

= (Z f (x)v *> Z ^ )v >) = Z Z = & M ^>- 

Definition 5.31. We say that A has a spectral gap a > iff 



//-• 



(4r,Aif/y H , for all ^ e I'd. 



(5.22) 



Theorem 5.32 (Spectral gap). A /zai a spectral gap a > if and only if (5.1) defines a bounded 



self-adjoint operator M : l L (X) — > wfr/z ||M|| < 



Proof. This follows immediately when either side of (5.21) is bounded from below by a > 0. □ 
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